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o ■ 
o : 

. Schiff moments were calculated for a set of nuclei with full account of core polar- 

Q I ization effects. A finite range P and T violating weak nucleon-nucleon interaction 

^ ' has been used in the calculations. While in the absence of core polarization the 

o : 

CN , Schiff moment depends on one combination of the weak interaction constants, in 

I the presence of core polarization the Schiff moment depends on all three constants 

^2 ' separately. The dominant contribution comes from isovector, AT = 1, part of the 

. weak interaction. The effects of core polarization were found to have in general a 

o : 

t;^ , large effect on the Schiff moments. 
O ; 

! PACS numbers: 21.10.Ky, 21.60.Iz 
C ■ I. INTRODUCTION 

■ - - ' time (T) reversal violating electro-magnetic (E-M) moments generated by nuclei. It is con- 
jectured that the dominant contribution to the atomic dipole moment comes from parity 
and time reversal violating components in the nucleon-nucleon (N-N) force. These compo- 
nents of the interaction induce P-T violating E-M moments in the nucleus. In turn these 
forbidden nuclear moments produce parity and time reversal violation in the E-M field which 
can induce static electric dipole moments in the atom. This way of reasoning is the working 
assumption in the study of the atomic dipole moment 

The measurements of the atomic dipole moment are presently very advanced and reached 
a very high degree of precision . These experiments are performed with neutral atoms. 
It has been shown years ago ^ that in such cases the field of the electric dipole moment of 
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the nucleus is screened by the field of the electrons and there is no effect when the atom is 
placed in an external electric field. This situation applies strictly only when one considers 
point-like particles. When the finite size of the nucleus is taken into account one can go 
to higher orders in the expansion of the E-M field and find that in the next order in the 
expansion it is the Schiff moment that can induce an electric dipole moment in the atom. 
(See next section). 

The upper limits for the existence of dipole moments obtained in the atomic measure- 
ments cannot provide directly limits on the presence of T and P violation in the nuclear 
Hamiltonian. One needs theoretical input. In order to make the connection one must have 
precise calculations of moments in nuclei produced by the T-P non-conserving part of the 
N-N interaction. 

In the past the theory of Schiff moments in spherical nuclei was limited to either a single- 
particle value (in Z-odd, N-even nuclei) or calculated for a simple two particle - one hole 
(2p-lh) configuration (in Z-even, N-odd nuclei). 

In this paper we will go beyond these simple calculations and include various effects of 
core polarization. As we will see these effects are large compared to the results obtained in 
lowest order. 

We should stress that an important aspect of this work is the use of a finite range 
interaction for the P-T-violating component of the N-N force. A finite range interaction 
was used recently in the calculation of the Schiff moment of ^^^Hg ^ and ^^^Ra[3]. In the 
present work we examine the influence of flnite range for a number of nuclei and discuss 
more general features that this leads to. 

We compute the Schiff moments for the following nuclei, for reasons explained below. 



A. Neutron-odd nuclei 

a. ^^^Hg: The most accurate upper limit on parity and time reversal violation has been 
obtained from the measurement of the atomic electric dipole moment in the Hg atom P|. 
The nucleus has (Z=80) has an odd number of neutrons. The lowest order contribution to 
the Schiff moment comes from a 2p-lh conflguration (a proton Ip-lh added to the neutron ) 
admixed into the ground state (g.s.). Because of the importance of this nucleus in the present 
experimental studies we have examined in more detail the influence of core polarization 
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effects on the Schiff moment . 

b. ^"^^Xe: The A=129 isotope of Xenon is also an even-Z, odd-N nucleus and the lowest 
order contribution to the Schiff moment comes from a Ip-lh proton excitation of the core. 
Similarly to ^^^Hg the atom of ^^^Xe is used in experimental studies of the atomic electric 
dipole moment ilj]. The upper limit for the dipole moment obtained in the experiment is 
only several times higher than in Hg. We have, therefore as in the case of Hg, found it useful 
to calculate more accurately the Schiff moment by introducing several effects coming from 
core polarization and not included in previous theoretical studies. 

c. "^^^Rn and "^^^Ra: The light isotopes of radon and radium are considered to be tran- 
sitional nuclei, still being close to spherical in their ground states 0,1^. The heavier isotopes 
of these elements (A=221-225) are found to have quadrupole plus octupole deformed shapes 
in the ground states 0, ^. It was suggested that in these heavier isotopes of Rn and Ra 
(as well as some other elements in this region ) the Schiff moments are strongly enhanced 
because of the quadrupole+octupole deformations 11|. It is therefore useful for the pur- 
pose of comparison to calculate accurately the Schiff moments in the neighboring, lighter 
isotopes where the onset of deformations has not yet occurred. 

d. "^"^^Ra: Here we calculate the Schiff moment assuming that the nucleus is spherical. 
We do this in order to have a reference value so that one can see the degree of enhancement 
when a full calculation , that includes deformations, is performed. 



B. Proton-odd nuclei 



In a Z-odd, N-even nucleus the odd proton can contribute directly to the Schiff moment. 
It is therefore usually assumed that the dominant contribution is the single-particle value of 
the last proton. 

e. ^^^Cs: We examined this point for ^^^Cs (Z=55) where the core-polarization correc- 
tion was evaluated together with the single-particle Schiff moment. 

f. "^"^^Fr: This nucleus (Z=87) was found to have an enhanced Schiff moment when the 

n 

quadrupole-|-octupole deformations were taken into account As in the case of ^^^Ra for 
the same reason we computed the single-particle and core polarization values taking '^'^^Fi 
to be spherical. 
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II. P AND T VIOLATING NUCLEON-NUCLEON INTERACTION 



The expectation value of the Schiff moment operator 



1 5 
S = — E e^l^^'r, - 3 < >ch ri) (1) 

can be non zero only in the presence of P and T violating nuclear forces. There are good 
reasons to assume that the exchange of a vr'^-meson is the most efficient mechanism of 
generating CP odd nuclear forces. This is due to the large value of the strong ttNN coupling 
constant Qs = 13.5 and to the small pion mass, as well as to the fact that outer proton and 
neutron orbitals in heavy nuclei are quite different. The P odd, T odd effective nNN 



Lagrangians are conveniently classified by their isotopic properties 

Ar = 0. Lq = go[V2{pn7r^ + npn') + {pp-nn)TT% (2) 

|Ar| = 1. Li = giNNTT^ = gi{pp +nn)n^; (3) 
|AT| = 2. L2 = g2[V2{pmv+ + np7T-) -2{pp-nn)'K°]. (4) 
The above Lagrangians generate a P and T odd nucleon-nucleon interaction of the form 

W{ra - Tb) = [(fi-oTa " Th + 92^7 a " Tfo - 3r„V^'^)) (cr„ - CT^) 

+gi{r:<Ta - r>b)] ■ V, , (5) 

rab 

where rrip is the proton mass and rab= | fa — | . 

Together with the finite range interaction Eq. (jSJ a phenomenological zero range effective 
interaction has often been used in estimates of the Schiff moment. It has the form 

G 1 

Wc{ra - n) = -^-^ iiVabCTa " VbaCTb) " Va5(r„ - r;,) + 



V'abi^a X a-h] ■ {(p„ - Pb), 5{ra - r^)}) , (6) 

where G is the Fermi constant. The interaction (jH)) can be obtained from (0) in the limit 
m-TT oo. In this limit we have the following correspondence between the interaction 
constants rjab and gi 

7gs{-go + 25-2 - 5-1), Vnn = 7r^5'.(-5'o + 25-2 + fl-l), 



Q^2js, . -^z ^x;, .,nn ^ . 
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Vnp = yr^9si9o-'2g2 + gi), Vpn = yr^9si9o - '^92 - 9i)- (7) 



^Js{9o-292 + 9i), Vpn = ^, 

The second term in Eq.® comes from the exchange matrix elements of the interaction (jS)) 
in the hmit 00 



\/2 \/2 

Vpp = Q^3s{90 + 2^2 + ^1), V'nn = Q^9si90 + 2^2 - ^l), 



V'pn = Vnp = TT-y 2^,(^0 - 92)- (8) 



Gml 

One should however remember that the correspondence given by (|7j) and (jH} is exact only in 
the limit of infinite pion mass. For finite pion mass this correspondence is somewhat condi- 
tional since the matrix elements of finite range interaction differ from the matrix elements 
of zero range interaction. Especially large is the difference between the exchange matrix 
elements. 

For a nucleus having an odd proton there are two kinds of contributions to the nuclear 
Schiff moment even in the absence of a strong residual interaction between the odd proton 
and the nucleons in the core One contribution comes from the odd proton. We call it 
a single particle contribution. Another contribution comes from the core nucleons. 

The weak interaction, Eq.(0), generates a weak correction to the nuclear mean field. The 
correction is 

,2 



SU{r) = ^{CT ■ n)T^ / r'W6(r,r')[(^7o-2^72)(Pp(r)-p„(r))- 
7rm„ Jo 



9iippir) + Pnir))], (9) 

where b{r, r') is a combination of spherical Bessel functions of an imaginary argument = 
^7r/2x/„+i/2(x), and kn{x) = ^j2/TixKn+i/2{x) 

b{r,r') = ii{mT^r)kQ{mT^r')6{r' — r) — iQ{mT^r')ki{mT^r)9{r — r'). 

Note, that the weak potential given by Eq.Q is proportional to r^. It is pure isovector. 
The contribution of AT = 0,2 channels to the weak potential is proportional to pp — pn, 
while the contribution of AT = 1 channel is proportional to Pp + Pn- It means that the 
weak potential is most sensitive to gi. The contributions of go and g2 are suppressed by the 
factor (A^ — Z)/A. This suppression has an important effect on the AT = channel. The 
finite range real A^A^-interaction ^ has both direct and exchange matrix elements. The 
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TABLE I: The contributions from the direct and exchange potentials to the single particle Schiff 
moment of ^^^Cs. The units are e fm^. 





9s9a 


9s9i 


9892 


direct 


0.011 


-0.109 


-0.022 


exchange 


-0.016 


-0.003 


-0.006 



direct matrix elements contribute to the potential in Eq.Q. For a finite range interaction 
the exchange matrix elements are usually small and in many cases can be safely omitted. 
For the AT = channel, due to suppression in the direct potential, the exchange and direct 
potentials become comparable. 

The weak potential Q produces a correction 5il)y{Y) to the single particle wave function 
ipu{r) of the odd proton. With this correction the expectation value of the Schiff moment is 

Ssp = {^Pu + S^,\S,\4j, + 5iP,) = {5^AS,\iP,) + (^.1^.15^.). (10) 

Note that this single particle contribution is absent for neutron odd nuclei. In Table 1 we 
show the single particle contributions to the Schiff moment of ^^^Cs from the direct and 
exchange potentials in all three channels. One can see that while in AT = 1 and 2 channels 
the exchange contribution is really small, in AT = channel it is even larger in absolute 
value than the zero order term and has opposite sign. The same behaviour was found also 
for 223pj^ nucleus. 

The weak potential 5U is defined by diagonal matrix elements of the two body weak 
interaction over core states. Another type of contribution to the Schiff moment comes 
from off diagonal matrix elements of the two body weak interaction. 

g (z/, 0|g.|lpl/i, lplh\W\Q, v) ^\W\\v\h, z/)(z/, lpl/i|g,|0, v) 

cove. / J j-p y— I • \ / 

Calculating the matrix elements in ()11|) we obtain the following equation for this contribution 

Score = Y.{^,^l\WW2,iy)^^^^^^^{'y2\S.Wl), (12) 



where are the occupation numbers and are the single particle energies. The states i^i and 
1^2 belong to protons, therefore, for odd protons nuclei Score is proportional to ripp while for 
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TABLE II: The single particle and the core contributions for two proton odd nuclei. Finite range 
weak interaction. The units are e fm^ 







9sgo 


9s9i 


9s92 


133Cs 




-0.005 


-0.112 


-0.028 




Score 


0.088 


0.088 


-0.176 


223pj. 


Ssp 


0.014 


0.084 


0.029 


223pj. 


Score 


-0.135 


-0.135 


0.27 



TABLE III: The single particle and the core contributions for two proton odd nuclei. Zero range 
weak interaction. The units are e fm^ 







9sgo 


9sgi 


9s92 


133Cs 


Ssp 


-0.06 


-0.23 


-0.0562 


133Cs 


^core 


0.18 


0.18 


-0.36 


223 


Ssp 


0.071 


0.19 


0.072 


223 JY 


Score 


-0.209 


-0.209 


0.418 



odd neutrons nuclei it is proportional to rjnp- The magnitudes of Ssp and Score are comparable 
for proton odd nuclei. This is in contrast to the P odd and T even nuclear anapole moment 
where the core contribution is smaller by a factor l/A^^^. In Table 2 we show the singe 
particle and the core contributions for two proton odd nuclei. The sum of three terms in 
the core contribution is proportional to the combination —go — gi + 2g2 ~ rjpp = —rjnp. 

At this place it is worth to compare the single-particle and the core contributions calcu- 
lated in the limit of zero range weak interaction. The results of calculations are presented in 
Table III. We see that the difference is approximatly by a factor 2. The matrix elements of 
a finite range interaction are smaller due to smaller overlap of the wave functions in matrix 
elements of the interaction. However, this is not a general rule. For the Hg nucleus, as 
we mentioned in the effect of a finite range is insignificant and goes even in opposite 
direction. For a finite range interaction we obtained 

Score = -0.085(5(s5'o + gsgi - '^QsQl), 
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while for zero range interaction 

Score = -0.058{gsgo + gsQi - 2gsg2)- 
For Rn nucleus the difference is 38% and for Xe nucleus it becomes again close to 2. 

III. CORE POLARIZATION 

The above results were obtained in a simple independent particle model (IPM). In our 
calculations we used a full single-particle spectrum including continuum. The single parti- 
cle basis was obtained using the partially self-consistent mean-field potential of The 
potential includes four terms. The isoscalar term is the standard Woods-Saxon potential 

= -l + exp[l-R)/ay ^^^^ 

with the parameters V = 52.03 MeV, R = 1.2709^^/^ ^^^^ ^ ^ o.742 fm. Two other 
terms Uis{r) and Ujir) were obtained in a self-consistent way using the two-body Migdal- 
type interaction |l7| for the spin-orbit and isovector part of the potential. The last term is 
the Coulomb potential of a uniformly charged sphere with Rc = I.ISA^^^ fm. 

In the next step, the residual strong interaction between the odd particle and the particles 
in the core should be taken into account. In IPM the state of an odd nucleus can be presented 
as 

k) = 4|o), 

where |0) is the ground state of an even core. When the interaction is switched on the odd 
nucleus state becomes a complicated superposition of the exited states of the core and the 
odd particle 

\u) = Aal\0) + B,,,{lplh)al,\lplh) + C,>,{2p2h)al,\2p2h) + ■■■. (14) 

(Note that the round brackets mean a perturbed state). This superposition can be conve- 
niently written as 

W) = Ual\0), (15) 

where If is a unitary transformation depending on a nucleon-nucleon interaction. The matrix 
elements of any operator O between the exact states of the odd nucleus can be expressed 
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by matrix elements of an effective operator O between the IPM states 

{v'\d\v) = {u'\d\u), (16) 

where O = WOU is the effective operator. For a one body operator, hke the Schiff moment, 
in first order in the nucleon-nucleon interaction the matrix elements of the effective operator 
coincide with Eq. ()12|l for the core contribution. 

^ rj — ri 

{u'\0\u) = {i^'\0\i^) + Y,{iy'i^i\VM _ ''^ V2\0\vi), (17) 

where V is the full nucleon-nucleon residual interaction. Leaving only Iplh intermediate 
states we can sum over perturbation series obtaining for the effective Schiff moment the 
equation 

{p'\S\u) = {v'\S\u) + Y.{p'vi\V\p2v) _ '^^ (z/2|S|z/i), (18) 

where (i/'|S|z/) is the single particle matrix element of the bare Schiff moment operator, 
Eq.(^. Eq .()18 p is just the RPA equation for the effective field in the terminology of Migdal's 
theory |l7|. This equation describes the effects of core polarization due to interaction of the 
odd particle with the particles in the core. 

Diagonal matrix elements of the Schiff moment are nonzero only in the presence of P and 
T violating interactions. For this reason the full nucleon nucleon interaction V must include 
this interaction 

r(ri,r2) = F(ri,r2)+Vr(ri-r2), 

where F is the strong residual interaction and W is given by Eq.Q. In the same way, the 
mean field potential is the sum of two terms U{r) + 6U{r), where U{r) is the main mean 
field created by the strong interaction and 6U (r) is the weak correction given by Eq. p()|l . As 
metioned above, 6U{r) creates the correction 6ip^{r) for every single particle wave function 
in all matrix elements in Eq. ()18|) . Since the corrections are small we can retain only 
the first order. To do this one should perform the substitutions V F + W and all 
"^u "^u + ^"^u in Eq. (jl8|) and keep the terms linear in W or 6tp. After this procedure the 
contributions to the Schiff moment can be written as a sum of three terms 



S = {6^,\S,m + i^ASzlS^u) + {iy\SS\u). 



(19) 
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The first two terms are those where the weak correction enters via the odd particle wave 
function. S satisfies the equation 

(zy'|S|z/) = (z/'|S|z/) + ^(z/Vi|F|z/2t/) _ "' HSlu,), (20) 

that differs from Eq. (jl8|) only in the interaction. Here only the residual strong interaction F 
enters in the equation. The effects of the weak interaction are entirely in the wave functions 
of the odd proton. Eq. ipUl) describes the well known effect of renormalization of nuclear 
moments due to coupling with particle-hole states in the core. Renormalization of the Schiff 



18(|. It was 



moment due to coupling with the isoscalar dipole modes were discussed in 
found there that the effect of renormalization for ^°^Bi nucleus is not significant, about 
15%. We found similar values both for Cs and Fr nuclei for renormalization of the isoscalar 
component of the Schiff moment. The isovector componet, however, undergoes significant 
renormalization. Due to strong repulsion in the isovector channel it leads to reduction of 
the isovector component. As a result, the single particle contrubution becomes reduced by 
the factor ~ 2. For Cs, instead of values cited in the first line in Table I we obtain for the 
renormalized single particle contribution Ssp 

Ssp = -O.OOS^fsS'o - 0.056gsgi - 0.016gsg2- 
The third term in (fT^ satisfies the equation 

77 — 77 

{u'\5S\u) = {u'\5SoW) + T.(^'''^\FM '^^ _ ''" {i^2\SS\u,), (21) 

that looks similar to ()20|) . However, the inhomogenious term (5So is completely different, 
namely 



UIV2 ^^1 

77 — 77 ~ ^ 

J2{u'u,\FM _ •'' m,,\S\u,) + (z/2|S|#.,)). (22) 

The first term on the rhs of Eq. ()22j) is the core contribution ()12|1 . where instead of the 
bare Schiff moment operator enters the renormalized operator S. The second and the 
third terms correspond to additional contributions where the weak interaction enters via 
corrections to the intermediate single particle states |z/i > and |z/2 >. Equations (20-22) 
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TABLE IV: Effects of core polarization on the Schiff moment for odd proton nuclei, for finite range 
weak interaction. The bare values of the Schiff moment, without core polarization, are listed in 
the first line for each nucleus. The units are e fm^ 





9s90 


gsQi 


9892 


133Cs 


0.08 


-0.02 


-0.21 


133Cs 


0.006 


-0.02 


-0.04 


223j;y 


-0.122 


-0.052 


0.300 


223 Pj. 


-0.009 


-0.016 


0.030 



describe all of the core polarization effects. The main contribution to (5So comes from the 
first term in Eq. This term depends on the renormalized Schiff moment S, therefore 

its contribution becomes reduced proportionally to the reduction of S, approximately by the 
factor ~ 2. The contribution of the second and the third terms in Eq.(j22I) is small. It is 
noticable only in T=l channel. 

In Tables IV, and V we present the results of calculations of the Schiff moment for sets of 
odd proton and odd neutron nuclei. The effects of the core polarization are twofold. First, 
without core polarization, for odd neutron nuclei the Schiff moment depended only on one 
constant rjnp (see Eq.((7j)). Now, the Schiff moment depends on all three constants. Second, 
the values of the Schiff moment are decreased. There are two sources of that reduction. 
First, it is mentioned above reduction in the first term in Eq. (l22j) . Further decrease comes 
from repulsive nature of spin-spin residual forces responsible for renormalization of (5So in 
Eq.dHD- 

IV. CONCLUSIONS 

We stress here two main conclusions from this work. First, the core-polarization con- 
tributions to the Schiff moments are sizeable. These corrections reduce the lowest order 
values obtained in the past by factors two to three. This has some important consequences 
when trying to interpret the implication of the experimental atomic dipole measurements. 
The reductions of the Schiff moments (compared to the values obtained in the past) means 
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TABLE V: Effects of core polarization on the Schiff moment for neutron odd nuclei, for finte range 
weak interaction. The bare values of the Schiff moment, without core polarization, are listed in 
the first line for each nucleus. The units are efm^. 





gsQQ 


9s9i 


9s92 


199jjg 


-0.09 


-0.09 


0.18 


199 jjg 


-0.00004 


-0.055 


0.009 


i29Xe 


0.06 


0.06 


-0.12 


i^^^Xe 


0.008 


0.006 


-0.009 


2iiRn 


-0.12 


-0.12 


0.24 


2iiRn 


-0.019 


0.061 


0.053 


2i3Ra 


-0.12 


-0.12 


0.24 


2i3Ra 


-0.012 


-0.021 


0.016 


225Ra 


0.08 


0.08 


-0.16 


225Ra 


0.033 


-0.037 


-0.046 



that using the experimental upper limits for the atomic dipole one obtains now a higher 
upper limit for the strength parameters in the weak parity and time reversal non- conserving 
component of the N-N interaction. Second, we have seen that the finite range reduces 
the mixing matrix elements between the opposite parity states in the case of spherical nu- 
clei. It was argued m that the large enhancement factors for the Schiff moments found 
in quadrupole-|-octupole nuclei [llj will be reduced when finite range interactions are used. 
We see that the enhancement factors (compared to the spherical case) will not be reduced 
because the effect of finite range has a similar effect in spherical and in deformed nuclei. 

We wish to thank V. Zelevinsky for helpful discussions. One of us (V.D.) wants to 
thank Tel- Aviv University for the hospitality. This work was supported by the US - Israel 
Binational Science Foundation. 
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